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ABSTRACT 

Aims. We studied the relationship between the average stellar abundance of several elements and the orbital evolution of stars in the neighbour- 
hood of the Sun. 

Methods. We used both observational data for 325 late-type dwarfs in a volume-complete sample and simulations of the orbital diffusion. 
Metallicities, ages, and initial position and velocities for the simulated stars are sampled from empirical distributions of these quantities in the 
Milky Way. 

Results. We found that that there is a relationship between the average stellar abundance of Fe, Na, Si, Ca, Ni, and Ba and the mean orbital 
radius of stars currently passing through the solar neighbourhood. The greater the difference between the mean orbital radius and the solar 
Galactocentric distance, the more deficient the star is, on average, in these chemical species. 

Conclusions. The stars that take a longer time to come from their birthplaces to arrive in the present solar neighbourhood are more likely to 
be more metal-poor than those that were born here. This result is a direct, independent indication that a tightly defined Galactic age-metallicity 
relation exists. 
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1. Introduction 



Although one of the core consequences of the chemical evolu- 
tion theory is the gradual increase in the metal content of the 
interstellar medium (ISM) and the progressive enrichment of 
subsequent stellar generations, some authours have found little, 
if any, indication that an age-metallicity relation (AM R) ex- 
ists amongst solar neighbourh ood late-type stars jFeltzing et alJ 
12001 1 Nordstrom et alJEo04h . None of these papers question 
the predictions of the chemical evolution theory concerning the 
enrichment of the ISM. The discrepancy between theory and 
data is usually interpreted as a result of an inefficient mix of 
the stellar ejecta into the surrounding medium, in spite of the 
growing evidence that the present dispersion in the ISM ove r 
distances of 1 kpc is very small JScalo & Elmegreerl 12004). 
However, other authors have found a much tighter relationship 
between age and metallicity ("see lRocha-Pinto et alJ200oT) . 

A way out of this deadlock was recently proposed by 
IPont & Everl f2004). who showed that the non-linear separation 
of the isochrones in the HR diagram can introduce a statistical 
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bias in the isochrone age measure ment, artificially as signing a 
greater isochrone age to the stars. IPont & Everl J2004I) showed 
that if a proper statistical method is used to correct for this bias, 
a tight AMR can be found from the same data used by some 
of the 'loose' AMR proponents. On the other hand, it is true 
that chromospheric ages also suffer from thei r own biases and 
cannot be considered perfect age estimators (Donahue, 1998; 
IPace & PaquinilEo04l) . Since much of the problem stems from 
the drawbacks of age determination methods, the existence of 
an AMR can be tested with the use of an age-related quantity 
that is independent of any age determination method. In this 
Letter, we show that the mean orbital radius of stars presently 
passing through the solar neighbourhood, R,„, is mildly related 
to the stellar age and that it can provide this independent test of 
the existence of the AMR. 



2. Sample selection and observations 

Our sample consists of stars from the Hipparcos catalogue 
( E^A , 1 19971) that also have uvby indice s measured bvlOlse 



1993i 1 19941) or compiled by lHauck &~ Mermillic 
( 1998) and which could be observed from the southern hemi- 
sphere. A heliocentric distance cut of 25 pc was used in order 
to have a sample of nearly 300 stars, so it would be possible to 
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complete the programme within a few observation nights at a 
small telescope. Our final sample has 325 stars. 

The observations were carried out at the Laboratorio 
Nacional de Astroffsica, as part of a survey of the chemical 
abundances of late-type dwarfs, during 15 nights between July 
1999 and October 2001. The coude spectrograph of the 1.60m 
telescope was used with the 1800 1/mm grating and a 1024- 
pixel CCD, pixel size being 24/y. The spectra were centred at 
/16145A, covered about 150 A, and had a nominal resolution 
of 0.3 A, with an average S/N ratio of 100 to 200. All spec- 
tra were reduced in the standard way; continuum normaliza- 
tion was performed homogeneously for all spectra. The spec- 
tral range contains a number of moderately strong, unblended 
transitions: 2 for Na I, 3 for Si I, 3 for Ca I, 10 for Fe I, 2 for Fe 
II, 5 for Ni I, and one for Ba II. The equivalent widths (EWs) of 
the metal lines were converted to elemental abunda nces using 
the MARCS model atmospheres described by E dvardsson e t al. 
l l 19931) . in a differential analysis with the Sun as standard star. 
Solar gf values were derived by forcing EWs measured off 
Moon spectra, obtained in the same conditions as the stellar 
spectra, to reproduce the standard solar abundances. Effective 
temperatures for the program stars were deri ved from (b — y), 
(B—V ) and (Bj- Vt) calibrations described by del Peloso et al. 
(2005). Surface gravities were obtained from the ionization 
equilibrium of the Fe I and Fe II lines. A full description of 
the abundance analysis will be presented elsewhere. The inter- 
nal uncertainty of the [X/H] abundances is ~ 0.1 dex, except 
for Ba II, which has more uncertain abundances, with only one 
spectral line measured. 

Radial velocities (v r ) for most of the 325 observed stars 
were calculated from Th-Ar comparison spectra. For a few 
of them, comparison spectra could not be obtained, and their 
v r were taken from the SIMBAD database. We calculated the 
U, V, and W components of the spatial velocity around the 
Galactic centre, from the v r , Hipparcos parallaxes and proper 
motions. Thereafter, the orbital parameters (mean orbital galac- 
tocentric radius, eccentricity, maximum heig ht above t he plane) 
for ea ch star were calculated as described in Rocha-Pinto et al. 
(2004, hereafter, R04). In this paper, we have adopted the 
present solar Galactocentric radius, R Q , as 8.0 kpc. 

A full chemokinematical analysis of this data sample will 
be presented elsewhere, but here, we are interested in an ap- 
parent trend between R m and the stellar abundances. Figure ^ 
shows plots of average elemental abundance, [X/H], averaged 
o ver several R m b ins, for Fe, Na, Si, Ca , Ni, an d Ba. According 
to lGrenonl ( fl987l) and lEdvardsson et alJ ll993). R„, can be taken 
as an indicator of the stellar Galactocentric birth radius, R*. 
IWielen et alJ Jl996h have shown that this relation between age 
and R ,„ can only be considered in a statistical sense. However, 
found a correlation between the age and R„, range for the 
stars presently crossing the solar neighbourhood: younger stel- 
lar generations have R m near R Q , while the older generations 
have the broader R m distributions typical of stars coming from 
different parts of the Galaxy into the solar neighbourhood. This 
finding is consistent with Grenon's proposition that R m stays 
more or less close to R* and reinforces previous claims that R m 
can be used to estimate R* . 
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Fig. 1. Average abundance of Fe, Na, Si, Ca, Ni, and Ba as 
a function of the mean orbital Galactocentric radius for stars 
presently crossing the solar neighbourhood. The average abun- 
dance is larger for stars having R,„ around the present solar 
Galactocentric radius, and it decreases by more than 0.2 dex 
for stars at extreme R m . The vertical and horizontal error bars 
are standard deviations around ([X/H]) and (R m ), respectively, 
at each bin. 

In Fig. n we show that the average abundances for stars 
presently crossing the solar neighbourhood also differ as a 
function of R m . The stars with R m w R Q have near-solar 
([X/H]), while stars supposedly born farther and farther from 
the solar neighbourhood (i.e., those with larger \R m - R Q \) have 
progressively decreasing ([X/H]). This decrease can be as high 
as 0.3 dex for Fe and 0.4 dex for Ba, between R m = 8 and 
R m = 6 kpc. 

An anonymous referee has pointe d out that the same trend 
can be seen in seen in Table 14 of lEdyardsson et al](ll993h : we 
have also found it for [Me/H] amongst late-type dwarfs from 
the Geneva-Copenhagen survey (see Fig.|2j. 

Disk stellar velocities peak near Cl(R*), the average rota- 
tional velocity at the stellar birth radius. Most of this motion 
goes in the azimuthal direction, not in the radial direction. 
Radial an d vertical mov ement is mostly driven by the orbital 
diffusion JWieletlfl977i) . on account of stellar encounters with 
massive objects like giant molecular clouds. As time goes on, 
the distance between the apogalactic and perigalactic radius in- 
creases due to the random kicks that stars get from collisions, 
and the chance the star could be caught at a different galacto- 
centric distance from its birth radius increases as well. Thus, 
since R m can be seen as a stellar birth radius indicator, we ex- 
pect that for our local sample \R m - R Q \ reflects the time the 
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star spent to come from R* to R G . In other words, \R m - R Q \ 
can be taken as an independent age estimator. It has an obvious 
advantage of neither depending on stellar evolution models nor 
being based on chromospheric properties. However, like other 
kinematical age estimators, it does not allow us to assign ages 
for individual stars, but to a stellar group. This is easily seen 
from Figure 7 by R04: The larger the age spectrum, and conse- 
quently the age uncertainty, the smaller \R,„ -R e \. Hence, while 
stars with R m « R have a wide age range, stars that came into 
the solar neighbourhood from other radii are, on average, older 
than the stars born near R Q . 

If this is the case, the relation between average abundance 
and R m seen in Fig. ^re fleets the existence of a relation between 
age and elemental abundances that is similar to the AMR; oth- 
er wise, if there were no AMR in the Galactic disk, as hinted 
bv lFeltzing et alJ l)200ll) . solar neighbourhood stars with large 
\Rm - R®\ (old stars according to the reasoning presented in 
the previous paragraph) could have any abundance within the 
typical abundance range of disk stars, and their average would 
not be dissimilar to the average abundance of stars having 
\R m - R e \ * 0. 

3. Simulation of the /^-abundance plane 

To test this hypothesis, we ran chemodynamical simulations of 
a stellar sample, using a reasonable Galactic model constrained 
by empirical data. Because we are interested in thin disk stars, 
we did not address the simulation of the thick disk or halo, al- 
though the effect of these Galactic components on the Galactic 
potential is taken into account. 

The starting position of a simulated star in X, Y, and Z coor- 
dinates with respect to the Galactic c entre is randomly ch osen 
following the thin disk density law of Sieg el et alJ (J2002). The 
star formation rate was considered to be constant, and a ran- 
domly uniform age from to 15 Gyr is assigned to each star. 
The simulated stellar metallicity follows the distribution 

[Fe/H] ~ n(z(t) + (R - R Q )d R [Fe/H],cr 2 ) , (1) 

where z(t) represents the AMR at age t, cr is the cosmic scatter, 
and <9«[Fe/H] the radial metallicity gradient at the galactocen- 
tric distance R. Finally, the starting velocities V r , V$, and V z for 
each star, in the radial, azimuthal, and vertical direction, respec- 
tively, with respect to the Galactic centre, were sampled from 
a normal distribution ~ NUVi), cr^\ where (V r ) = (V z ) = 
and (K)) = ®(R ), i.e., the rotation velocity at R, taken from 
Law et alJ J2005h . We took [o>,., cr v o>J = [23, 12,9] km/s 
fromlR04[ 

The orbit of each star was then followed from its birth — 
given by the randomly assigned age — to the present time , 
using a program developed by Chris Flynn (details in R04). 
To simulate stellar encounters with massive objects, a random 
push in each velocity component was given at each integration 
step. The magnitude of the push in each velocity component 
was set by requiring that the increase in the velocity dispersion 
with age of the stellar group for the simulated stars in the so- 
lar vicinity fitted the observed age-velocity dispersion relation 
(given by R04). After integrating the orbit, stars whose present- 
time galactocentric radius differed by less than 0.5 kpc from 
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Fig. 2. Comparison between the observed and predicted R m - 
[Fe/H] trend for stars passing through the solar neighbourhood. 
The hatched grey area and white squares correspond to stan- 
dard deviations around the average [Fe/H] for the same R m in- 
tervals for 13017 late-type dw arfs in the Geneva-Copenhagen 
survey JNordstrom et alll2004l) . The lines correspond to the re- 
sults of a Galactic model with no AMR (dashed line), loose 
AMR (dotted line), and tight AMR (solid line), as detailed by 
the parametrizations in Tabled 

Table 1. Parametrizations used in the simulations 



function 


parametrization 


no AMR 


z(t) = -0.15; cr = 0.3 dex 


loose AMR 


z(t) = 0.0239 - 0.0355t - 0.0007r 2 




cr = 0.22 dex 


tight AMR 


z(t) = 1.51 -44.25(34.23 -r)- 1 




cr = 0.1 dex 


single-slope 8r [Fe/H] 
two-slope flflfFe/H] 


-0.07 dex kpc" 1 

-0.05 dex kpc 1 , R < R Q 

-0.03 dex kpc -1 , R > R e 



R Q were selected as "solar neighbourhood" stars. It is from this 
simulated subsample of stars passing through the solar neigh- 
bourhood that we take the predicted 7?,„-abundance relation to 
compare with the results from Fig.^ 

Since our goal was simply to show that the observed trend 
in the /? m -abundance plane is determined by the existence of a 
well-defined AMR, we used three single parametrizations for 
the AMR: i) pure scatter relation, for which no change in the 
average metallicity or abundance scatter occurs with time; ii) 
a loose AMR, correspon ding to a polinomial fit to the data in 
lEdvardsson et alJ Jl993l) : iii) a relatively tight relation accord- 
ing to which the average metallicity of younger stars is higher 
than that for older stars. The adopted parametrizations for these 
two cases are given in Tabled 

Our simulations consider an initial number of 6000 stars, 
in order to have nearly 300 stars within 0.5 kpc from R G at the 
end of the integration to compare with the observational sam- 
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pie. Fig. |2 summarises our simulations of the /?,„-abundance 
plane. Three curves are shown, which correspond to the R m - 
<[Fe/H]> trend found using the tight AMR, loose AMR, and 
no-AMR case, each with the two-slope metallicity gradient 
given in Tabled Clearly the no-AMR simulation failed to pro- 
duce the arc-shaped R m -([Fe/H]) trend seen in the observa- 
tional data. This failure is not difficult to understand, because 
if there is no AMR in the Galactic disk, there is no relation be- 
tween how much time a star takes to come from R* to R G and 
its metallicity. The R m -([Fe/H]) trend for this case is close to 
the adopted average metallicity for this parametrization ([Fe/H] 
as -0.15 dex). The trend also decreases very slowly with R m on 
account of the adopted 3s[Fe/H]. 

The solid curve in Fig. |3| is very similar to the observed 
trend. It is the resulting R m -([Fe/H]) trend for the case of a 
tight AMR with very small intrinsic metallicity dispersion. We 
have made no attempt to find a best-fit curve, since the goal 
of the simulations is to show that the observational trend be- 
tween R„, and [Fe/H] can be explained with reasonable as- 
sumptions for the disk parameters, provided that there is a tight 
AMR. Several other parametrizations for the AMR and metal- 
licity gradient were used, in order to verify the effect of the 
parametrization on the R m -([Fe /H]) trend. We verified that the 
shape of the resulting curve is most sensitive to the overall 
metallicity growth from the beginning of the Galactic disk to 
the present time, as well as to cr. A lower overall metallicity 
growth and/or larger cr flatten the simulated R m -{ [Fe/H]) rela- 
tion, making it look closer to the no-AMR case. On the other 
hand, the effect of <9ft[Fe/H] is more pronounced on the tails 
of the j?„,-([Fe /H]) curve. We ver ified that a two-slope metal- 
licity gradient ( M aciel et all 120061) yields a better eye-ball fit 
to the observed trend than a single-slope gradient. However, 
the effect is small and should be investigated more before any 
conclusions can be made. 

Our simulations have only considered the R„-([Fe/H]) 
trend because there is much more information available in the 
literature for the observed relations between [Fe/H], t, and R. 
However, the same exercise could be done for Ca, Na, Si, Ni, 
and Ba, as long as their abundances with respect to H have 
grown since the beginning of the Galactic disk and the intrinsic 
cosmic scatter at any age was low. It is interesting to verify in 
Fig.[Qthat the biggest difference in the abundance at R m « R Q 
and large \R m -R e \ is found for Ba. Barium is thought to present 
a larg er overall growth with age than iron (Edvard sson et all 
1 1993b . since the bulk of its nuclesynthesis is due to long-lived 
AGB stars, which make their presence felt in the chemical evo- 
lution of the Galaxy with a lag with respect to iron sources. 
From our simulations, a larger growth in abundance with age 
increases the arc-like shape of the R m -{ [Fe/H]) relation, on ac- 
count of the difference in abundance from the older stars with 
larger \R„, -R e \ to the average abundance stars with R m « R G . It 
is also interesting to note that Ca and Si (and perhaps also Na) 
seem to have the flatter profiles with R m in Fig. 1 . That they are 
a-elements, almost exclusively products of type II supernovae, 
could account for this, since a less well-defined age-abundance 
relation is expected in this case. 



4. Conclusions 

Spectroscopic abundances for 6 elements in 325 stars show 
that, on average, solar neighbourhood stars with large \R m - R Q \ 
have lower abundances, regardless of the chemical element 
considered. This relation is interpreted as evidence of a well- 
defined age-abundance relation, since stars with large \R m —R Q \ 
would mainly be old. A number of simulations of the orbits of 
disk stars born with chemokinematical properties that are rea- 
sonable for the Galactic disk confirm this interpretation. 

The constraint provided by our data on the tightness of the 
AMR is relatively strong. Figure [5] shows that a loose AMR, 
with cr as 0.20 dex, is still compatible with an arc-shaped R m - 
([Fe/H]), albeit much less pronounced than what is observed. 

Our method highlights the possibility of studying the AMR 
and its dispersion through a technique that does not depend on 
age determination drawbacks. This could be achieved through 
a larger stellar sample and accurate (preferably, spectroscopic) 
abundances in order to better define the empirical R„,-{ [X/H]) 
relation. 
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